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This study proposed an enhancement technique for improvising the 
estimation technique in iterative back projection (IBP) by using the 
Lorentzian error function with a sharp infinite symmetrical filter (SISEF). 
The IBP estimation is an iteratively based error correction that can minimize 
the error reconstruction significantly. However, the IBP has a drawback in 
that it suffers from jaggy and ringing artifacts as a result of the iterative 
reconstruction method and the absence of edge guidance. Furthermore, 
because the IBP estimator tended to oscillate at the same solution frequently, 
numerous iterations were required. Therefore, this study proposed edge 
enhancement to enhance the estimator by using the combination of the IBP 
with Lorentzian SISEF to produce a finer high-resolution output image. As a 
result, the SISEF is used to improvise the estimator by providing high 
accuracy of edge detail information for enhancing the edge image. At the 


same time, the Lorentzian error norm helps to increase the robustness of the 
IBP algorithm from contamination of additional noise and the ringing 
artifacts. 
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1. INTRODUCTION 

Most digital imaging applications demand higher-quality image acquisition. The finer quality of the 
image or high resolution (HR) image is crucial to gain the most useful information as much possible [1]-[17]. 
However, in the reality of image acquisition system was introduced a degraded quality of the output image or 
low resolution (LR) image by the hardware limitation. This low-quality image was degraded by many factors 
contributed including internal and external environments. For instance, a finite number of resolution 
acquisition sensors may draw limited information on the digital image. Moreover, the motion blur appeared 
on the image when the poor handling of the acquisition device between the scenes. Besides that, the 
acquisition system may introduce a noisy image due to limitations of the sensor’s device and additional 
signals during the transmission medium [18]—[33]. However, the hardware improvement not always reliable 
in most cases. It because the improvement of the sensor acquisition system would lead to a large expenditure 
budget, and unreliable for real applications. Thus, image recovery through the super resolution technique is 
the best choice to produce a finer quality image with enlarge the image size capability. 

This study is intended to highlight the constraints on the super resolution (SR) technique and 
provide some improvements to improve the SR. The essential in SR is the reconstruction technique, which is 
based on the manipulation of the global imaging problem model. This model is generated based on the 


Journal homepage: http://ijece.iaescore.com 


2540 O ISSN: 2088-8708 


existence of several limitation factors in the real acquisition system such as constraints that occurred in the 
point spread function (PSF) of the camera that contains information of blur, translation, and decimation 
factors [34]. In addition, some constraints on the image acquisition may lead to the presence of global noise 
that contaminated the acquired image [21]—[23], [26], [29]-[33]. The SR reconstruction is referring to the 
recovery techniques in image processing that obtains a high resolution (HR) image from observed single or 
multiple LR images. To produce the HR image, the researcher has manipulated the mathematical expression 
of the imaging model which describes the sensor characteristics, such as motion blur, decimation, and shifted 
translation in the camera PSF [1], [2], [10], [11], [15]. 

The SR reconstruction method consists of three basic tasks: registration, interpolation, and 
restoration [35]. This task can be implemented separately or simultaneously which depends on how the 
technique is adopted. The re-construction process takes sub-pixel shifted information from the sequence of 
LR images or sensor parameters provided by the image registration task. Indeed, the accuracy of sub-pixel 
motion is a very important factor for SR reconstruction solution to produce a finer image [21], [23], [27], 
[31], [34], [36]-[47]. Then, the interpolation task is taking place to align the nonuniform space of the LR 
image onto a uniformly HR image grid. Finally, the reconstruction task deploys to restore the finer HR image 
by removing the blur and noise effects. 

The iterative back projection (IBP) is among the first spatial-based SR reconstruction. This 
technique is formulated by Irani and Peleg [41], [46], [48], this approach is similar to the reconstruction of 
2-D objects from 1-D projections in computer-aided tomography (CAT). Each of the LR image resolution 
pixels is a projection of a region in the scene whose size is determined by the imaging blur. Then high 
resolution is determined using a similar back-projection method in CAT. The technique is equipped with sub- 
pixel accuracy as image registration to align the LR images on the HR grid before increases the size image 
resolution. Then, the HR image is estimated by back projecting the error between simulated LR images and 
observed LR images. The gradient or error is determined the same as the least square error in l, norm 
computation. Moreover, the IBP technique requires an initial guess for very first-timer estimation. If the 
initial guess estimated is close to the optimal point in the search space, then the estimation will converge 
smoothly. 

This back-projection technique can minimize the error of reconstruction efficiently in an iterative 
manner. To correct the next estimated image, this error reconstruction is back-projected to the HR grid. 
However, by removing the blur effect, this technique can improve the image intensity [24], [33], [49]-[52]. 
The iterative procedure, on the other hand, contains defects that cause ringing artifacts at the edges. This 
condition emerges as a result of incorrect back projecting due to the lack of edge guidance during projection 
[22], [33], [49], [53]-[57]. Therefore, this study proposed the improvement with the edge enhancement 
technique embedded with the IBP reconstruction technique by using a combination of the Lorentzian 
function and sharp infinite symmetrical exponential filter (LSISEF). The SISEF function has high 
localization edge detection and it is used to preserve the possible edge detected. Coincide, supplies the edge 
guidance during the back projecting to avoid across edge projection, thus the ringing effects be reduced. 
Besides that, the Lorentzian error norm is deployed to suppress the amount of error in the cost function. As a 
result, the suggested methodology can improve HR estimation by maintaining lost frequency information 
while being resistant to ringing effects and extra noise [32], [58]—[60]. 


2. SUPER RESOLUTION 

Since the conventional super resolution (SR) approach suffers from HR estimate imperfections. This 
study gets a lot of attention from the researcher to produce a better technique for improving the ordinary one. 
Therefore, Chiang and Boult [61] proposed that the resampling algorithm be improved by incorporating the 
image warping methodology into the IBP reconstruction. It was accomplished by combining the imaging 
system’s degradation model with integrating resampled data, which resulted in a more accurate estimation of 
IBP. Cohen and Dinstein [62] indicated that the final image resolution might be improved by using a 
polyphase filter to replace the geometric transformation inside the IBP algorithm. Furthermore, Kim et al. 
[63] suggested an edge enhancement strategy that uses anisotropic diffusion to minimize noise and highlight 
edges. To remove blur efficiently, Bose et al. [64] proposed optimal regularization using L-Curve for 
handling the degraded image. To provide a better reconstruction technique, Rajan and Chaudhuri [65] used 
generalized interpolation scheme development for image resizing and proposed a combination magnetic 
resonance fingerprinting (MRF) and MAP estimator for SR reconstruction. Sun ef al. [66] improved the 
quality hallucinate image by implementing Markov-chain based inference algorithm inside back-projection 
Rajan and Chaudhuri [67] proposed super-resolution reconstruction using an machine learning (ML) 
estimator with MRF regularization to estimate the depth map and focus on the scene image. 
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Furthermore, the resolution improvement via implementing the 3D motion is estimated recursively 
in an extended Kalman filter and by a novel image warping procedure [68]. The IBP reconstruction error can 
be reduced by incorporating adaptively nonlocal back projection [69]. To improve the reconstruction errors 
reconstruction in the IBP, Dong et al. [70] proposed the non-local IBP algorithm in which the enlargement 
process incorporates adaptively the non-local information. The resolution improvement in efficient highlight 
edges image can be done by embedding the un-sharpening mask by using the SRUM algorithm [71]. 
Furthermore, Yan et al. [72] proposed a combination Papoulis-Gerchberg extrapolation inside the IBP to 
enhance the spiral CT slice image. Furthermore, the IBP technique's efficiency is improved by using a 
wavelet locally adaptive approach to estimate the first initial guess rather than the traditional interpolation 
approach [73]. Following that, Liang and Gan [74] show that using the non-local IBP fast technique, they can 
improve edge identification during the initial interpolated image. Besides that, the L*a*b domain was 
proposed by Bengtsson et al. [75] for improving the human visual system in the HR image. By maintaining 
the edge information details, the output HR image will be improved. 

Besides that, to recover a high-resolution image, Bareja and Modi [76] proposed combining the IBP 
technique with the canny edge detector to reconstruct a high-resolution image. As a result, the error 
difference between degraded and estimated images is used to recover lost frequency. Moreover, the proposed 
hybrid technique was introduced between the IBP technique with the Cuckoo optimizer by improving the 
initial guess and HR output image [77]. Rasti et al. [78] proposed an approach via using the interpolation 
method on the LR before registered inside the IBP algorithm. Cheref and Yousfi [25] used the anisotropic 
diffusion approach to preserve possible edge pictures. In SR reconstruction, Maiseli et al. [79] used adaptive 
Perona Malik as edge preservation. 

This work continued with improving the IBP result by reducing the ringing effects by applied 
variant and local variance variant based on adaptive analysis [80]. The visual quality image improved while 
applied the sparse linear regression and iterative back projection. Additionally, they modified the Gaussian 
high pass filter to refining the initial guess [81]. When the residual error was reduced by utilizing anisotropic 
diffusion to assist the IBP estimator in producing a sharper HR image with edge preservation, the accuracy of 
the reconstruction was improve. Hassen and Jahmeerbacus [82] proposed combining the IBP with a Canny 
edge detector and a Gabor filter to improve back-projection error minimization. Furthermore, Nayak and 
Patra [83] proposed using the P-Spline and MuCSO-QPSO algorithms to improve the edge image of the IBP 
methodology, with the P-Spline providing crisp edge information and MuCSO-QPSO optimizing the 
estimation approach. 

In work [84], a simultaneous task in SR reconstruction is introduced for robust outliers and 
preserving the edge image by implementing the Huber norm for robust regularization. Meanwhile, in [83], 
[85] the proposed a robust iterative super resolution reconstruction by using Lorentzian Thikonov with a 
Bayesian MAP estimator. This technique is used to remove multi kinds of artifacts and improved the 
estimation [83], [85]. The IBP algorithm also suffers from ringing effects. This problem can be reduced after 
employing total variation as a penalty function during the reconstruction process [86]. The study for 
removing the noise and ringing effects is continued in work [87], through employing an effective edge- 
guided interpolation method with a bilateral filter for significantly removed the artifacts. The regularization 
method proposed by Feng and Lei [88] is a combination of the total generalized variation (TGV) and Shearlet 
Transform to eliminate the artifacts and preserve the edges of the image. To suppress the ringing artifacts in 
the IBP output, Yang et al. [89] introduced a variation and local variance variation based on an adaptive 
inside the reconstruction. Many proposed algorithms have shown improvements in the robustness towards 
outliers and noise. However, the improvements should bring benefits to the SR quality image and the 
computational complexity as well. 


3. METHOD 
3.1. Imaging model 

This section is discussed the PSF model of the SR reconstruction method. The digital application 
represents an image in matrix form as a column-wise lexicographically order matrix notation. Which the 
input is given a number N observed of LR images {Y;,}/-,, where each of the input images has matrix size 
[M,, X Mg]. Then the reconstruction output HR image X has matrix size [L x L], where commonly the size 
L > M, and L = æM, the a is a multiplication enlargement factor. Since each of the observed input image 
has arbitrary geometric warping, linear space-variant blurring, and uniform rational decimating done on the 
ideal HR image X. Furthermore, the input images assumed each of them are degraded by non-homogeneous 
additive noise. Thus, all parameters involved above description to an analytical model show as in (1): 


Yk = DkBkWkX + Ny (1) 
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where, Wp is [L? x L?] the matrix representing geometric warping performed on the image X. Then, B,is the 
linear space-variant blur matrix of size [L? x L?] and Dx is a [M? x L?] matrix size represents as a decimation 
operator. Where 7, stands for the additive noise in the kth measurement with the positive definite 
autocorrelation matrix with matrix size [M, X Mg]. During the experiment setup, all matrices (Dk, By, Ty, Nx) 
are assumed to be known in advance after preparing by the prerequisite image registration stage. 

The parameter matrix Dx represents the decimation factor, which this ratio factor represents an 
enlargement factor size between the ideal reference image and the kth measurement image. This ratio is 
obtained from the division between the size of a pixel in the measured image Mẹ” and the reference image 
L?. If choosing a very high ratio factor of the decimation and will cause the ill-posed problem and the image 
suffers from lost a lot of information. 


3.2. Iterative back projection 

The IBP technique's fundamental idea is to reconstruct the HR image from several degraded input 
images that have been detected. The difference between simulated and observed LR pictures is projected 
back to estimate the difference . The estimating method is repeated until the HR picture generates the same 
LR images as observed LR images after passing through the same blurring and downsampling method. The 
IBP method can divide into two phases are simulates the observed image and back projection the error to 
reconstruct the HR estimation [46]. The first phase is the algorithm has to construct the simulated LR 
images Y,. 


Yk = D,By WX (2) 


By referring to (2), the observed LR image known as Y with the number kt” of images. While the production 
of pixel intensity of LR images is influenced by the degradation function derived from the HR image X. This 
degradation function kernel is the yield a combination of the blur kernel function Bx and translation wrapping 
function W;, before downsampling by the decimation operator D. 

Initially, the IBP technique procedures starting with deploys an average method to compute the 
initial guess X HR image. The next stage is the imaging process produces the set of simulated LR images 
which are known as yo However, initial guessing may result in multiple solutions rather than a single 
correct solution. The algorithm's ability to discover a unique, quicker, and smooth answer is dependent on the 
accuracy of the first estimate technique. Recently, the average of the observed registered LR images, as 
shown in (3), is a good initial guess method [46]. 


XO = E Diaa We p) (3) 


When the initial guess HR image X© is equal with the reference image, thus the simulated images 
Y¥,° = DyB,W;.X should be matching with the observed images Y; and no further correction is needed. If 
the between Y, and Y;°still not identical then, by back projecting the error and combining the current 
information in the HR grid, the initial guess is improved even further. This step of the process is yielding a 
better result to the next HR image X © such as in (4): 


tH = fh + yal Qr [Yk — Qf] (4) 
where 
Qk = Dy By We (5) 


This procedure keeps repeating iteratively until the error is minimized enough. The error function can be 
described as (6). 


e= Puun- ni © 


3.3. Infinite symmetrical exponential filter (isef) 

The Shen Castan edge detector or infinite symmetrical exponential filter (ISEF) is an ideal 
smoothing filter with a relatively simple recursive method [90]. The implementation of the differentiation 
optimal exponential in discrete form is the heart of this filter. It has a positive impact on identifying the edge 
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picture since it responds less to noise and has a high degree of precision in edge localization. Edge detection 
was proposed using the maxima of gradient (GEF) or the zeros crossing of the second directional derivative 
along the gradient direction (SDEF). When it comes to edge detection, the Laplacian Gaussian filter 
confronts a tradeoff between noise insensibility and precise localization. 

To overcome this challenge, the ISEF created a linear filter based on a one-step model and 
multi-edge detection. This optimal smoothing filter is an asymmetric exponential filter with a very simple 
recursive algorithm for infinitely large window sizes. The difference between the input and output of this 
recursive filter, which is capable of detecting edges with low noise sensitivity and high localization precision 
[90]. For symmetric exponential filters, the ISEF filter generates the GEF and SDEF from the first and 
second derivatives of the operator. The (7) is the normalized symmetric exponential filter in one dimension: 


fi) = CA) + fp) — aô (x)) (7) 
where 


C= —— (8) 


(2-ao) 
The exponential is the asymmetrical function which can be distinguished into the right side function 


f(x) and the left side function f)(x). Thus, the second derivative of exponential function for the SDEF 
along the x-axis can describe in (9): 


ð? 
Lay) = zall y) * fy) (9) 
where, in the discrete form of the SDEF shown in (10): 


Ll2, y) = y) * iO”) * fy) * (A) O) -2U y) * AAO) * fy) (10) 


The equation above also can be duplicated for calculating the gradient along the y axis direction. For 
instance, this equation flows can be illustrated in the process flow in Figure 1. 


1, (x,y) I.2(%Y) 


Figure 1. ISEF process flow 


3.4. Lorentzian function 
The noise outliers can reduce by deploying a robust error norm in the estimation stage. Then, the 
Lorentzian error norm is deployed to handle this constraint. This error norm which detected by Lorentzian 
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then projecting back those errors between an estimated HR image with observed LR images on the next HR 
estimation process [91]-[93]. Based on this characteristic, thus the Patanavijit proposed the estimation SR 
reconstruction using combination with the Lorentzian error norm. Where the Prop is the Lorentzian error 
norm function and this parameter value compute in pixel-wise operation such in (11): 


pron) = log h+] (11) 


where, the T is the Lorentzian constant parameter also known as the soft threshold value. When the value 
x < T then the function obeys the norm L2 form, but the result be saturated when inverse condition, where 
Wror is known as the Lorentzian norm influence function. The Lorentzian error norm function and 
Lorentzian influence function are illustrated in Figure 2(a) and Figure 2(b) respectively. 


2x 


Pior) = saa (12) 


reer 
2TG +x 


PLOR (Xx) 


(b) 


Figure 2. Lorentzian response (a) Lorentzian error norm and (b) influence function [83] 


3.5. Propose ibp technique with lorentzian sharp isef regulator 

The performance of the edge enhancement technique depends on the accuracy of the edge detector. 
The finest edge detector could offer high precision localization, a small error rate, and only detect the real 
edges. This study realizes the ISEF gives a promising accuracy during edge detection. This filter uses an 
exponential function based on designing the filter. Whereas this filter provides better edge detection than the 
filter Gaussian based. This study prefers to use the second derivative filter SDEF in the ISEF technique. It 
uses the zero-crossing approach during obtaining edge detection. The SDEF provides the edge guided with 
higher precise edge localization [90]. From the observation, the edge detected by SDEF has a wider contour. 
To increase the precision localization and decrease the error detection rate. This study has proposed sharper 
edge detection by subtracting the different edges of initial guess with the edge of current estimation such as 
in (13). 


Vf = (f — ft) * ISEF — (fe) * ISEF (13) 


This output of this difference provides a thinner ad sharper edge detection. As a result, this estimator 
could determine the HR much accurately and faster than previous methods such result in the next section. 
Furthermore, to increase the robustness and computational efficiency, this technique collaborates with the 
Lorentzian norm, such as shown in (14). 


fest = fe +Z Piort" + TWior Cf — fe) * ISEF — (fi) * ISEF] (14) 


The Lorentzian norm function Pior is used to improve the error residual computation and the 
Lorentzian influence function Yor increases the gradient detection value. While the threshold level T in the 
Lorentzian norm since T > 1 and it allows rapid correction in the residual error of the estimation. This 
condition catalyzes the estimator to reach the local minimum faster. Besides that, while the threshold level of 
the Lorentzian influence function T, > 1, the level of gradient detected was increasing, and the estimator 
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capable to preserve many possible of the edge’s information for each iteration. This proposed technique’s 
process flow can be illustrated in Figure 3. 


oo 


Set minimum error and maximum iteration 


r"> min Error 
Or t< max iteration 


Update HR image 


Calculate LR simulated gf = Hf” 
Compute residual error rë 


ISEF Filter 


Estimate the solution 


1 
fy =f" +57 Pror(r") + t(rorl(f — fr) * ISEF — (fe) * ISEF]) 


Figure 3. The IBP with Lorentzian ISEF regulator process flow 


4. RESULT 

The results of monitoring are shown in Tables 1 and 2, respectively. This section presents the 
discussion on experimental results have obtained by several techniques has explained before and including 
the proposed IBP with edge enhancement technique. The MATLAB platform is used as an experiment 
simulator in these experiments. Aside from that, these tests employ a regular image from the ISI and 
DESCAI databases, with a resolution of 256x256 pixels. 

The initial procedure of this experiment creates a sequence of four LR images obtained from the 
image database. Each of the LR images created must have shifted from the original image by a pixel in the 
vertical direction. After that, they have applied the effect of camera PSF. The effect involves the process of 
blurry and decimation operation. This shifted LR image was then convoluted with a 3x3 Gaussian low pass 
filter with a standard deviation of 3. In each direction, the image was down sampled by a factor of two. As a 
result, four LR images are created using the same technique with varying motion vectors in vertical and 
horizontal dimensions. The PSNR magnitude is used to evaluate the suggested technique's performance. In 
addition, the cost computation complexity indicator is determined by the number of iterations required to 
complete the HR calculation. 
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4.1. Situation 1: The performance of the proposed technique with the variation of the input image 

Table 1 reveals that the proposed LSISEF regulator methodology is responsible for the majority of 
the PSNR results. This means the proposed technique provides a better edge enhancement technique 
compared to other techniques. By looking at the average results, the LSISEF provides an improvement of 
about 2.12 dB from the conventional IBP technique. The key to improvement depends on the feature of the 
edge detector that provides a high localization precision, high rate detection, and small error rate detection. 
This positive change is illustrated in Figure 4 that provided a comparison output with the proposed technique 
which Figure 4(a) illustrates an input LR image and Figure 4(b) represents output HR produced by IBP and 
Figure 4(c) shows the proposed IBP with LSISEF technique output. Meanwhile, Figure 5 illustrates the 
overall output HR image from various techniques, such as Figure 5(a) represents the output of combination 
IBP with Canny method and Figure 5(b) shows the output by using patel method. In addition, Figure 5(c) 
represents the output of the cheref method and Figure 5(d) illustrates the output of the proposed technique 
which combined the IBP and LSISEF technique. 


Table 1. Measurement of PSNR for IBP edge enhancement 


IBP The IBP technique with 
Canny Patel Cheref LSISEF 
LENA 27.730 28.018 28.651 28.034 28.990 
CAMERAMAN 26.477 26.758 27.220 27.319 27.321 
MANDRILL 23.284 23.336 23.464 23.486 23.494 
PEPPER 25.941 26.391 27.956 29.023 30.676 
CLOWN 24.560 25.067 26.673 27.932 29.868 
JETPLANE 26.020 26.353 27.177 27.571 27.707 
BRAIN 28.041 28.699 30.131 31.319 32.681 
SHIP 28.923 29.188 29.494 29.714 29.609 
USAF 22.670 22.656 22.183 22.624 22.783 
BRIDGE 23.982 24.178 24.648 24.824 25.763 
Average 25.763 26.064 26.760 27.185 27.889 


LSISEF 


(a) (b) (c) 


Figure 4. Comparison output with the proposed technique (a) input LR image, (b) HR image of IBP 
technique, and (c) HR image LSISEF (proposed) 


CANNY CHEREF LSISEF 


(a) (b) 


(d) 


Figure 5. Overall output HR image (a) IBP+canny, (b) patel method, (c) cheref method, and (d) IBP+LSISEF 
(proposed) 


The LSISEF edge detector has higher localization precision than the other edge detector, therefore 
in most cases of the image, the proposed ISEF technique gives better performance. The Lorentzian norm 
function and influence function can increase the level of edge detection, thus improves the correction task 
during estimation. The saturation level in the Lorentzian function is set equal to 1.5, consequence the norm 
function and influence function of the Lorentzian act as a weight of residual error and gradient detected, 
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respectively. This process increases the amount of gradient and relights the edge of the image. As a result, 
while our proposed methodology is capable of preserving high information detail losses and producing a 
finer output image with good performance, it has difficulties handling crucial images such as those with a lot 
of high-frequency information. The Mandrill image, for example, has a lot of high-frequency data, such as 
hair surrounding the Mandrill face. When the LSISEF detector simply traced the image's edges in the 
horizontal and vertical directions, this disadvantage arose. As a result, this detector failed to offer details of 
high-frequency information in some directions. 

On the other hand, in some circumstances, Cheref’s methodology also delivers the highest quality. 
This approach delivers the best quality, with a gain of roughly 1.42 decibels. This method worked effectively 
when dealing with challenging photos like Mandrill and USAF pictures. Cheref's anisotropic diffusion can 
offer edge detail in a variety of orientations (north, south, west, east). As a result, this method can provide 
detailed edge information as well as a promising result in a challenging image. Furthermore, this edge 
detector has a low sensitivity to reduction edges and has a good rate of detection. Thus, this technique is 
among the best edge enhancement technique. However, this technique has less localization compared to the 
LSISEF technique. 

Shreyas [94] have proposed the combination of the IBP technique with the ISEF edge detector. As 
aforementioned before, the ISEF provides high precision in edge localization and less response to the error. 
Therefore, it can increase the quality image up to 0.997 dB in average value. Unfortunately, this edge 
detection only coverage about vertical and horizontal directions, because of this limitation the improvement 
only preserves the edge in limits directions. At the same time, the edge detected by this technique is not tiny 
enough and it may invite small error detection. Consequently, the Patel works failed to provide a higher 
improvement on the IBP reconstruction technique. 

The Canny method is Gaussian based technique, and the result showed this technique only provides 
just a little improvement [76]. This technique has less error detection since it provides a very sharp edge 
image detected. However, this filter has disadvantages in the sensitivity detection rate. When a Gaussian- 
based filter has a contradiction function between giving a smooth function and an accurate localization rate, 
this problem occurs. As result this technique only capable to increase about 0.301 dB and the lowest 
improvement compared with other edge enhancement techniques. 


4.2. Situation 2: Blurry effect information unknown 

The objectives of this experiment are to evaluate the robustness of the IBP edge enhancement 
technique toward the blind blurry effects. This experiment assumes the image registration stage is failed to 
predict blur effect information accurately from the degraded input image. To imitate the real situation, this 
experiment uses injecting the unknown level of Gaussian blur effect to multiple input LR images. The blurry 
effect has 5x5 matrix size is used with varying values of sigma from 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0. To 
perform the de-blurring process, the kernel system is required to pick a random sigma value. Thus, this 
experiment uses the small sigma value is 0.5 as the kernel system random default sigma value. This 
experiment measuring the quality and cost computation of the IBP with the edge enhancement. 

Generally, the collection of the result in Table 2 shows that the proposed techniques capable to 
manage the blind blurring condition. Specifically, this technique successfully recovers the lost information 
due to the effects of low to high blur levels applied to the input image. The larger standard deviation value of 
blur effects may cause a lot of high-frequency information to become attenuated and lost. In addition, the 
previous IBP method has lacked determination on preserving the losses. Thus, the IBP failed to estimate the 
HR image accurately. The result shows the IBP method capable to provide the highest quality output if the 
blurry information used is the same as a blurry amount in the degraded input image. However, the 
performance of IBP going to drop if blurry information applied diverts far from an actual blurry amount in 
the LR image. For example, in the case of the standard deviation of Gaussian blur equal 0.5 is applied to the 
input image and the IBP kernel system. As a result, this technique produced a quality image of 33.594 dB. 
However, the quality of the output image decreased to 22.187 dB when increases the amount of blurring up 
to 3 which is the standard deviation value that has deviated far from the system expected. 

This IBP weakness can improve by implementing the edge enhancement technique. The result in 
Table 2 shows, the LSISEF technique produces the highest quality images compared to other techniques. In 
average result, this combination gives improvement up to 3.182 dB compared to conventional IBP. Besides 
that, the average result of the SISEF method is capable to increase the quality image by about 3.073 dB better 
than the ordinary IBP technique. This result also conveyed the ISEF based is capable to provide good edge 
detection in highly blurry situations. In addition, the sharping technique on the proposed technique also 
contributed to provide good results. However, the quality performance is slightly decreased from the Cheref 
method. This condition happens, when the clipper function tends to remove some amount of high-frequency 
information in the Anisotropic Diffusion. As a result, this technique slightly is degraded the performance 
quality of the Cheref method by about 0.539 dB. 
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However, the edge detector of Gaussian filter-based techniques such as Canny regulators failed to 
preserve the highly blurred edge image. This technique has a less rate of sensitive detection of the possible 
edge compared to anisotropic diffusion and ISEF. Therefore, during the sigma blur=3.0, most of the 
techniques have a weakness to predict the HR image and difficulty in recovering the lost information at the 
same time. Besides that, the ISEF-based edge enhancement method provides a solution to solve this 
difficulty. Therefore, the proposed technique and technique proposed by Patel give a promising result during 
handling the effect of highly blurring. 

In addition, this technique successfully provides the highest quality PSNR whenever the other 
technique failed to do so. This experiment proves the LSISEF can resist the condition of the highly blurred 
image since this image lost much high-frequency information. However, the edge detector of the Gaussian 
filter-based technique failed to preserve the weakest edge image. Therefore, most of them not efficient during 
handling the highest blur image. Moreover, during the standard deviation blur is 3, most of the technique has 
a disadvantage to predicting the HR image and coincide difficulty in recovering the lost information. 
However, the LSISEF based edge enhancement method has a solution to solve this difficulty. Therefore, the 
proposed technique and technique proposed by Patel give a trusty result during the highly blurry effect given. 
In the future, it is planning to utilize this LSISEF method to improve the effectiveness of quality threshold 
ARTMAP [95], [96] in pattern recognition problems. In addition, further investigation will be done to 
increase the performance of the moment invariant technique [97] for the feature extraction process with the 
use of the LSISEF technique 


Table 2. The performance comparison of the combination IBP with edge enhancement toward the blind 


blurry effect 
IBP The IBP technique with 
Canny Patel Cheref LSISEF 
Sigma = 0.5 33.594 33.592 32.528 33.151 33.046 
Sigma = 1.0 27.436 27.436 28.613 27.400 28.711 
Sigma = 1.5 24.611 24.611 27.919 25.281 27.982 
Sigma =2.0 23.24 23.240 27.699 24.415 27.751 
Sigma = 2.5 22.561 22.561 27.600 24.009 27.648 
Sigma = 3.0 22.187 22.187 27.548 23.789 27.594 
Average 25.60483 25.6045 28.65117 26.34083 28.78867 


5. CONCLUSION 

This study proposes using the Lorentzian Sharp ISEF (LSISEF) regularization as a robust edge 
enhancement to improve the IBP reconstruction technique problem. In the past, the IBP reconstruction 
methodology was plagued with ringing artifacts, necessitating a high number of iterations to complete the 
estimation. Furthermore, due to fuzzy effects, this approach loses the potential to preserve lost 
high-frequency information. As a result, in this study, the Lorentzian Sharp ISEF regulator was designed to 
provide a remedy to the IBP reconstruction defects by providing precise edge enhancement with light cost 
computation. As a result, this proposed technique provided 2.12 dB better PSNR compared to the 
conventional IBP technique after examined with various input images. Additionally, this proposed technique 
provided 3.182 dB PSNR better than conventional IBP in a vary of blurry PSF functions. This result came 
about as a result of the method's high precision localization and low error rate detection when used as an edge 
detector. These benefits enable the estimator to calculate the corrective process more precisely and quickly. 
At the same time, the Lorentzian norm and influence function increase the quality residual computation and 
gradient detection, respectively. 
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